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ABSTRACT 


IMPLEMENTATION OF AND MEASUREMENT WITH 
THE LIPA-TECHNIQUE IN A SUBSONIC JET 


LIPA (Laser Induced Photochemical Anemometry) was used to measure velocity, 
vorticity, Reynolds stress and turbulent intensity distributions in a subsonic jet. The jet 
region of interest was the area close to the jet-orifice. The LIPA-technique is a non- 
intrusive quantitative flow visualization technique, consisting of tracking a phosphor- 
escing grid of fluid particles, which is impressed by laser-beams directed into the flow. 
The phosphorescence of biacetyl gas was used to enable tracking of the impressed light 
grid. In order to perform measurements in a jet, LIPA was developed and implemented 
for the specific flow requirements. Nitrogen was used as the carrier gas to avoid 
quenching of the phosphorescent radiation of the tracer gas biacetyl by ambient oxygen. 
The use of sulfur dioxide to sensitize phosphorescent emission of biacetyl was 
examined. Preliminary data was used in a discussion of the potential of the LIPA- 
technique. 
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LIST OF SYMBOLS 


Roman 
A Area 
Acj Biacetyl 

C Circumference 
CpQ Proportionality constants 

CR Confidence region 
E Energy 

E, Triplet energy 

FDC Fluid dynamic quantaties 
f Focus length 
I Intensity 

I e Intensity of emission 

If Intensity of fluorescence emission 

Ip Intensity of phosphorescence emission 

Ij Intensity of laser line energy 

k Rate constant 

l d Diffusion length 

m Grid mesh size 

p Pressure 

r Jet exit nozzle radius 

RCR Relative confidence region 
mas Root mean square 

RSE Relative standart error 

S Path along which circulation is calculated (circumference of a grid box) 
SE Standard error 

s Distance of movement of grid intersection points 

t t-factor (see error analysis) 

u Fluctuating part of streamwise local velocity 

u’ Root mean square of u 

U Local mean velocity 

U c Core velocity 

v Fluctuating part of normal local velocity 

v* Root mean square of v 

x Streamwise Cartesian coordinate 
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x, Coordinate along laserbeam 

y Normal Cartesian coordinate 

Gre ek 

a Absorption coefficient 

r Circulation 

e Extinction coefficient 

O p Phosphorescent quantum yield (emitted intensity/absorbed intensity) 

At Time interval between laser pulse and picture detection 

At, Time interval between laser trigger out TTL- signal and actual laser pulse 

x e Emission lifetime 

z { Fluorescent lifetime 

x p Phosphorescent lifetime 

to z Vorticity component perpendicular to grid plane 

Observed emission length in streamwise direction (still camera experiments) 

Superscripts 

Indicates relation to 'reference grid' 

Indicates relation to 'distorted grid' 

1 Indicates singlet state 

3 Indicates triplet state 

Other 

[ ] Indicates mol-density (mol/volume) 

Mean value of local quantity 
< > Mean value of local quantity (same as overbar) 
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1. Introduction 

Optical techniques are widely used in fluid mechanics to observe and measure 
properties of flow fields such as velocities and densities. Many of these techniques are 
qualitative but of great value in guiding intuition for further research by quantitative 
means. Optical techniques are usually known for their largely nonintrusive properties as 
compared with methods like the Pitot tube or the hot-wire technique. 

Flow- visualization techniques may be coarsely subdivided into two categories: 
those that make use of light scattered by tiny particles in the fluid and those that make 
use of variation in fluid properties (i. e. refractive index). Among the methods that rely 
on scattered light, laser Doppler anemometry is now a standard means of obtaining fluid 
velocities. In laser Doppler anemometry, the fluid velocity can be measured with high 
accuracy as a function of time but only at a single point in the fluid at any given time. 
Particle Image Velocimetry (PIV) has been developed to enable researchers to look at 
instantaneous spatial data. The ultimate aim, of course, is the simultaneous determination 
of fluid velocities in a whole volume of a fluid. First steps in this direction have been 
taken with the development of speckle photography, which can give the instantaneous 
velocity field over a complete plane of interest in a fluid. Among the methods that make 
use of refractive index variation, holographic methods are established. All of the above 
techniques require the use of particles in the flow. 

LIP A ( Laser Induced Photochemical Anemometry) has been recently developed 
and has the ability of simultaneous determination of fluid dynamic quantities like veloc- 
ity over a whole volume of a fluid. This technique makes use of excited luminescent 
particles or molecules, thus, special seed molecules have to be introduced into the flow as 
the normal lifetime of the excited states is not long enough with respect to the time scale 
of the flow. Such molecules are the key component of the LIPA measuring technique. 

LIPA can also be performed for liquid flows using photochromic dyes or 
luminescent particles, which are premixed with the fluid and excited by a laser at specific 
locations in the test section. Due to the low vapor pressure of photochromic dyes, and the 
static charge that builds up in luminescence particles, these laser-marking methods are 
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not well suited for gas flows. One way to mark fluid elements in gases in a similar fash- 
ion as done with liquids would be to use laser-induced phosphorescence. In contrast with 
fluorescence, which involves quantum-allowed transitions and short lifetimes (10- 6 -10- 

8 s), phosphorescence occurs via forbidden transitions and, therefore, has a long radiative 
lifetime (>l(Hs). 

LIPA was previously conducted successfully in media such as water and kerosene 
[1,2] where luminescent chemicals can be dissolved easily. Tests in a Mach 3 jet gas flow 
[3] showed the potential for the LIPA investigation of a jet flow, however, no calibrated 
data was obtained. The jet was seeded with droplets of a water diluted luminescent 
chemical. The first successful LIPA investigation of a gas flow field has been made re- 
cently in a motored two-stroke engine during scavenging [4]. 

Biacetyl vapor, the luminescent chemical that is used in the present investigation, 
was first introduced as a flow field diagnostic material by Epstein [5] who made use of 
the fluorescence of biacetyl. Whereas biacetyls fluorescence has been used later on by 
McKenzie [6], Hiller [7] and Hilbert [4] applied biacetyls phosphorescence to flow field 
diagnostic. Because biacetyls phosphorescence, in contrast to its fluorescence, is strongly 

quenched by oxygen, the latter implementations have been made with nitrogen-biacetyl 
gas mixtures. 

There have been a large number of investigations of turbulent jets mainly because 
of their engineering importance, simplicity and relevance to other turbulent flows. Two 
books have been published which are largely devoted to this subject (see references [8] 
and [9]). Most experimental velocity measurements in turbulent jets have been made with 
hot-wire instrumentation. However, other studies have recently been performed using 
laser-Doppler velocimetry (LDV). Studies in supersonic jets are more difficult, because 
probes interfere with the shock structure. The particles needed for LDA may also cause 
problems. Therefore, LIPA, as developed here, may one day be a powerful tool for 
supersonic jet research. 

A general jet flow may be separated into three regions, each having a distinctly 
different character. One region consists of a potential core of uniform mean velocity and 
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low intensity of turbulent fluctuations 
bounded by a shear layer. This region starts 
at the nozzle exit and has a length of 4 or 5 
nozzle diameters. The mixing region, which 
lies between the central core and the 
undisturbed gas of the surroundings, where 
the velocity gradients are large and the 
intensity of turbulence is high, is a second 
region. The third region, where the central 
core and the mixing region blend into a 
completely turbulent gasstream starts at about 
8 nozzle diameter behind the nozzle exit. 


In this work LIPA is applied in an axial symmetrical round jet of a nitrogen-bi- 
acetyl gas mixture. The diameter of the jet exit nozzle is 10mm and the Reynolds number 
(based on the inner orifice diameter) of the room temperature (20°C) jet flow is about 
12400. The 47 measuring points employed here (see figure 1) cover an area of a approxi- 
mate distance from the jet centerline ±3r and a distance from the jet nozzle in streamwise 
direction of 0 to 7r. According to the estimated boundary of the potential core with Har- 
shas correlation (see [10], pg.: 20) and of the mixing zone with the potential theory, 23 
measuring points lie in gas of the surroundings, 14 measuring points are in the mixing 
area and 10 measuring points are in the area of the potential core. 
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distance from jet nozzle, x/r 

Figure 1: Position of the 47 measuring 
points 
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2. Experimental set-np 

A typical LIPA experimental set-up consists of three main components. One is 
the flow creation, which includes the seeding, if necessary, of the carrier gas with the 
luminescent chemical and the devices to generate the flow field of interest. A second 
element creates the light beam grid and consists of a pulsed light source (usually a pulsed 
laser), mirrors, lenses and a beam splitter to direct the light to the measurement area and 
beam dividers to create a grid pattern. A third component, that detects and stores the grid 
patterns (commonly a movie, CCD or CID camera with controllable shutter timing, 
which stores the detected information on a media like regular film material, video tape or 

disk), has to be, either case, coupled with the light-source to control the delay time be- 
tween light pulse and detection. 


The chosen set-up for the LIPA-implication for measurements in a free jet is 
schematically shown in figure 4. 

The light source is a Lambda Physiks LPX 220 pulsed excimer laser. It is 
charged with XeCl gas and the buffer-gas is neon. It emits ultra violet light at 308 nm. 
The initial beam size is about 5 mm by 20 mm and the maximum pulse energy that can 
be emitted during a pulse-length of 20 ns is 220 mJ. The maximum pulse energy de- 
creases with increase of the pulse repetition rate, that can be controlled internally (up to 
100 pulses/second) or externally (up to 250 pulses/second). The laser beam is directed by 
broad band aluminum coated mirrors and a 50:50 308nm dielectrically coated beam 
splitter through bi-convex quartz lenses of a focal length of 300mm to the beam dividers 
in order to create a grid of laser beams in the jet area of interest. The average pulse 
energy carried by one grid beam is determined to be about 3 % of the initial pulse 

energy, whereas a grid beam diameter (without scattering by the flow) of about 0.5 mm 
is detected. 

The implementation of the LIPA-technique requires 'lines’ of laser light to create 
a grid pattern. A beam divider is used to create several beams from one initial beam. In 
order to lose as little light as possible a beam divider that is shown schematically in fig- 
ure 2 is employed. The design of the beam divider resembles an oversized blazed reflec- 
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tion grating. Apart from diffraction 
and absorption effects, it provides no 
loss of incident laser energy if the 
angle between the incident beam and 
the metal base ( 9 ) is adjusted properly. 

Since the applied broad-band excimer 
laser radiation shows only a low degree 
of coherence, additional diffraction 
patterns (which could be observed if 
HeNe laser light was used for the 
irradiation) are absent so that the grid 
of beams in the test area is well defined. The aluminium 
coating of the mirrors is optimized for shallow incidence 
angles. Mesh sizes of 4.2 mm to 5.2 mm are generated in 
the actual experiment. 


The flowfield is generated by a 17m/s nitrogen jet 
(nozzle exit diameter: 1 cm) seeded with approximately 5 
mol% biacetyl, issuing into a 12 cm diameter and 74 cm 
long plexiglass tube. A slow side jet of the same gas 
mixture is provided to prevent self-feeding of the main jet. 

Quartz plexiglass windows are fitted into the tube to let the 
laser beams pass through and a honeycomb in the tube exit 
create a pressure gradient in order to reduce incoming oxygen and hence quenching of 
biacetyl emission by oxygen. The tank, that is connected to the exit nozzle, is provided 
with honeycombs (see appendix F) in order to reduce swirls in the initial jet. Seeding of 
the jet fluid with biacetyl (liquid at room temperature and atmosphere pressure) is 
accomplished by room temperature gas-dispersion (see figure 3). Nitrogen of 99.995% 
purity, provided by an industrial gas tank and pressure reduced by an industrial regulator, 
is dispersed through liquid biacetyl in a mixing chamber to evaporate biacetyl. To 
disperse the nitrogen through the biacetyl, a open ended tube is chosen in order to reduce 
the pressure in the nitrogen delivery system [ 1 1 ]. 
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Figure 3: Gas-dispers- 
ion bottle 
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Figure 2: Beam divider to create parallel 
beams (specular reflection) from one inci- 
dent beam 
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A gated ICCD (intensified charge coupled device) camera (ICCD-576 manufac- 
tured by Princeton Instruments) is applied to detect the grid patterns. Its detector 
(Thomson-CSF TH7883FO-2: 576 columns x 384 rows; spectral range of 400-1060nm 
sensitivity; dynamic range of 14-bits) is operated by a detector controller (Princeton 
Instruments ST- 130) that interfaces to a personal computer (Dell 325D). The camera is 
equipped with a 58mm lens (f-stop: 1:1.2) and a 52mm close-up lens. This optical set-up 
provides for the picture area a resolution of 0.0735 mm/pixel. The excimer-laser is the 
source of the system trigger. Two pulse generators are employed to set the delay between 
the laser pulse and the picture detection. Pulse generator #1 (Princeton Instruments FG- 
100) controls the time delay between the excimer-laser trigger signal and the actual laser 
pulse (At,=300ns), the time duration of the gating and prevents gating through CCD 

readout to reduce background noise. Pulse generator #2 provides pulse generator #1 with 
the actual delay between laser pulse and detector gating (At), an additional delay on At,. 

A second pulse generator is necessaiy because of the limited delay range of pulse 
generator #1 (20-1700 ns, whereas a delay time of At=30 to 90ps is required). A two 
channel oscilloscope is used to display and measure the time delays. 

Two types of datagnds are taken; undistorted and distorted grids. The undistorted 
grids are captured when the laser is firing, whereas the distorted grids are captured at a 
specified time delay after the laser pulse. Because of the low image detection repetition 
rate of the camera assembly (about 3 pictures/minutes), the distorted and undistorted 
grids are detected in different events (not during the same illumination process). Since 
indeed the undistorted grid pattern is detected at the time of the laser pulse, it is created 
by straight lines, and hence fundamentally a reference grid for all the distorted grids that 
are taken with the same experimental set-up. 

The camera assembly stores the raw data (pictures) on the hard drive of the per- 
sonal computer. In order to download the data sets into a Megavision 1024 XM image 
processor the data has to be reformatted to a 8-bit 512x512 pixel format (512 of the 
original 576 columns, 8 of the original 14 bits). This is necessary, because the image 
processor can only read data of a certain format. 
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Figure 4: LIPA-experiment set-up for measurments i: 
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3. Experimental p rocedure 


3x 1 The LfPA techniqnp. 


3x 1,1 Reyiny of the measurement te/'finiqnr 

The technique of Laser Induced Photochemical Anemometry (LIP A) makes use 
of ex, ted photochromic chemicals to measure important fluid quantities such as velocity 
and vorticity over a two-dimensional area of a fluid flow. 


Excited photochromic molecules can emit light for a certain time-interval (x e ) af- 
ter excitation. Hence, a photochromic molecules containing flow-volume radiated by a 
laser pulse (pulse length « Tj can emit light for a certain time-interval after the laser 
pulse. If non intrusive, with the flow moving photochromic chemicals are used, the path 

of a marked flow volume can be followed during the lifetime of the impressed light 
emission. 


The LIPA technique marks and follows points in the flow-field by employing 
photochromic molecules. The flow is seeded with a photochromic chemical and lines of 
pulsed laser light are directed into die flow. This impresses a grid of intersecting, light 
emitting fluid lines within die flow. The lines and especially die intersections of the lines, 
which define points in the flow, can be followed for a certain time after its creation. 

Two successive pictures of the grid are required to obtain fluid dynamic quanti- 
ties of the flow field. Separate pictures taken at die time of die laser pulse (a undefomied 
re erence grid') and after a time delay shorter than the lifetime of the emission (a 
deformed gnd'). By measuring the distance and direction each intersection travels and by 
knowing the time delay between each photograph, die two velocity components in the 
gn plane, turbulence intensities, Reynolds stress and vorticity can be calculated. 
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3.1,2 Algorithms tP obtain fluid dyn amic quantities 

To describe the procedure to obtain fluid dynamic quantities from a 'reference 

grid' and a 'distorted grid’, one 'grid box' is 
chosen (see figure 5a). If the mesh size of the 
grid is small enough, one box can be thought 
of as a fluid 'particle'. 


A fluid particle moving in a general 
three-dimensional flow field may have mo- 
tions about all three coordinate axes. The 
LIPA technique, in means of the application as 
it is shown here, is limited to the projection of 
a three-dimensional motion of the particle onto 
the plane of the photos which is parallel to the 
initial plane of the grid. 

As a particle moves with the flow, it 
may undergo several motions that can trans- 
late, rotate and deform the 'grid box' in the 
plane of the photos. Because the history of this 
specific fluid particle is known, the displace- 
ment Asj = (x i -x i .,y r -y i .) ^ the velocity 
Uj = As j / At of each comer i can be calcu- 
lated. The latter are average velocities over the 
time interval At, thus, they are designated to 
an average gnd box* (see figure 5 b and figure 
6: (1,2, 3, 4)) whose comers are positioned at 
the midpoint between the associated comers of the reference and distorted 'grid box' 
( x i’X) = ((Xi-+V)/2,(y i .+y i .)/2). 



figure 5: Procedure for conversa- 
tion of the displacement between 
the reference (l’,2',3’,4') and the 
| distorted (1",2",3",4") 'grid box’ to 
velocities. 
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To calculate the vorticity from the 
'grid box' (1,2,3,4) and the velocities (Ui) 
the definition of the circulation T is used: 



Figure 6: Decomposition of the velo- 
cities used to calculate the circulation 
around a 'grid box' of the area A and the 
circumference C. 


r = |0»dS (3.1) 

c 

To relate the surface integral to an 
area integral, Gauss’ theorem is employed, 
and leads to: 

r = J • n dA (3.2) 

A 

<?v 

with: (O z = — — (3.3) 

ax ay 


The velocity vector at each comer is converted to the directions of the box side- 
lines (see figure 6). Assuming a linear change of the velocity components along the box 
sidelines, U • dS in equation (3.1) is estimated for every box side by forming the average 
of the comer velocity components in the box sideline direction: 


r = |u«ds= 

C 


Ju«dS+]u»dS+Ju*dS+Ju*dS 

12 3 4 


(3.4) 


£ Uu±0^ . s u + 0,1 + 03,, . ^ + * 0,3 . + °4 - l ± ♦ §4., 


with: 



The average vorticity component at the centroid of the fluid particle normal to the 
picture plane cd z can be obtained by dividing the calculated circulation T by the area A of 

the fluid particle (equation 3.2): 
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= X (3-5) 

This technique has the advantages of avoiding a second differencing of the 
experiment data. 


By following this procedure for the other mesh elements, the vorticity compo- 
nents normal to the picture plane, and though in the plane of the initial grid, can be ob- 
tained at many locations over a two-dimensional field in a flow simultaneously. 
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3.1,3 Error analysis 

An error analysis, which indicates the expected relative error associated with the 
fluid dynamic quantities is shown. This error analysis is based on the assumption, that 
errors are generated only by inaccurate evaluation of the position of the grid intersection 
points. Other possible errors origins are neglected in this analysis. 

Assume that there is a uncertainty of the evaluated position of an intersection 
point (reference grid: 8s'; distorted grid: 8s") that moves a distance As during a time 
interval At. The expected errors for the single local measurement (single intersection 
point of single grid) are calculated as [12,13]: 

local velocity: 


£U X = <5U y = Su = Sv - 


j(Ssf+(Ss'f 

At 


(3.6) 


local Reynolds stress: 


8- (uv) _ ICAflEt - If***) 2 | 
uv llvuj IvJ fluj V v J 


(3.7) 


= <5Ur *. 






At 


iVuJ VvJ 


local turbulent intensity: 


u 2 / Uc^ = <?(u 2 ) 2 ^ Z ^ SU X j(&f +(&") 2 t 2 


u 2 /U c 2 


u 


At 


u 


(3.8) 


The uncertainty in the measured local vorticity is approximated by introducing 
the average, vorticity generating, velocity Un* along the sidelines of a gridbox. This 
velocity can be calculated by the average of the differences of the values of the two 
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velocities along opposide grid sidelines divided by 2(see also chapter 3.1.2). Because this 

velocity is the average of four independent velocities, its error and the error in the mesh 
size m is estimated by: 


J 4 *f t*M’) 2+ (&") 2 T 

L 

At 


(3.9) 


1 

= — * 
2 


f(Ssf + (Ss"f 
At 


1 fin 
2*ypz* At 


Assuming a square gridbox of the area A=m 2 , the relative error in the vorticity is 
estimated to (see also equation 3.4): 


local vorticity: 



= ^(&’) 2 +(<5s") 2 * 



The error analysis shows that the fluid dynamic quantities that can be calculated 
from the location of the single intersection points, involve different uncertainties. 
Whereas the uncertainty of the local velocities is only a function of the accuracy in find- 
ing the intersection points of the undistorted and distorted grid and of the time delay At, 
the uncertainty of Reynolds stress and turbulent intensity is also a function of the 
expected local fluctuation in the velocity. 
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The uncertainty in the local vorticity is approximated by a function of the mesh 
size, a average velocity along the sidelines of the grid box, the reading error and the time 
delay At. A classical error analysis that indicates the expected error resulting from double 
differencing a grid to obtain vorticity is presented in reference [14]. To compare the 
performed relative error estimation of equation 3.10 to the classical analysis, equal 
assumptions are made. During a time interval At one sideline of a square grid of the 
meshsize m is displaced by m/2 in sideline direction. This leads to U S j de =m/(8*At) and 
© z =l/(2*At) in equation 3.10. 


Denoting the reading error 


^(i5s') 2 + (<5s ") 2 as 5s and replacing in equation 3.10 


yields to: 



(3.10a) 


The result of the analysis in reference [14] is: 


S(t) z = 



(3.10b) 


A higher constant is obtained with the present error analysis, however, both 
equations contain the same function of At, 5s and m. 

Since cameras of all kind are basically limited in their resolution (instrumental 
limitation), a minimum uncertainty of the measurement is always involved. The maxi- 
mum accuracy intersection points can be determined is half of the resolution. The here 
employed camera detector provides a maximum resolution of 384 x 576 pixel which re- 
solves the comer velocities of a 10% deformed, imaginable 384 x 384 pixel mesh, by 2% 
maximum. 

Improve in the measurement accuracy of the velocity and hence in the other fluid 
mechanic properties can generally be achieved by increasing the time delay At, if the 
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resolution is unchanged. However, this requires larger mesh sizes and involves less local 
and temporal resolution of the fluid dynamic quantities. Increase in the resolution of the 
detector may decrease the minimum uncertainty but is only expedient if the actual read- 
ing accuracy of the intersection points is expected to be more precise than the current 
detector resolution permits. 

The reliability of the value that is the average of fluid dynamic quantities calcu- 
lated from single measurements is certainly larger, if mainly random errors occur (here: 
determination of intersection points by hand). Denoting the maximum occurring error in 
the reading of the intersection points by 5s’ m „ and 6s" m „, the standard deviation is esti- 
mated by 5s ,-2/3 * 5s mM and 5s ,=2/3 * 5s" m)lx [13]. If n is the number of measure- 
ments (distorted grids), then the (relative) standard error (R)SE of the fluid dynamic 
quantities is: 


standard error of local velocity: 


SE(Ux) = <5Ux = ~== * /( &, ») 2 + (&"«)* 


At 


2 


3*/tTT 


At 


relative standard error of local Reynolds stress: 

2 )/(<%> max) +(& >l max) 2 


(3.11) 


(3.12) 


3*^n^T 


At 
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relative standard error of local turbulent intensity! 



(3.13) 


relative standard error of local vorticity: 


RSE(fii z ) = — i 

COz 
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Because the expected local fluctuation of the velocities does not change with the 
local averaging of the data, the standard error in the mean of all the shown fluid dynamic 
quantities alters with the same function of the number of measurements. Note u' and v’ 
instead of u and v in the error equations of the mean local Reynolds stress and turbulent 
intensity. The root mean square quantities have been introduced because the mean of the 
fluctuation part of the velocities is zero by definition. 


The confidence region (a quantity for the mean reliability) of the mean value is a 
function of the standard error, of t (a factor that is dependent on the required statistical 
safety and the number of measurements; see [15,16] and appendix E for t-distribution) 
and of the number of measurements. The deviation between the mean and the real value 
decreases with the number of measurements and decreases the confidence region of the 
mean value. The relative confidence region (RCR; confidence region: CR) of the fluid 

dynamic quantities (FDQ) as a function of the relative standard error (RSE; standard 
error: SE): 


(R)CR fdq — ~^== * (R)SEfdq 


(3.15) 
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To achieve a statistical reliability of 95% that the real fluid dynamic quantity is in 
the confidence region of the mean of 24 measurements of a fluid dynamic quantity 
(t=2.08, like for the performed experiments), the (relative) confidence region is: 

(R)CR fdq = 0.43 * (R)SEfdq (3.16) 

Measuring over an area of large mean velocity gradients may involve large rela- 
tive uncertainty as well as large relative confidence region gradients over the measuring 
area if a constant time delay At is used. Assuming a constant uncertainty of evaluating 
the grid intersection points and a constant relative velocity fluctuation in the hole area, 
the absolute velocity fluctuations are small in areas of low mean velocity. Therefore, the 
relative error and the relative confidence region of the calculated fluid dynamic quantities 
is, respectively, proportional to 1/As, hence to the reciprocal of the local mean velocity. 
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3.2 Implementation of the LIPA-technique in a free iet 

3.2.1 Tracer gas, carrier gas and stimul ation process selection 

Because the turbulent length scales are small, a gaseous tracer is required for the 
implementation of the LIPA-technique in a free jet to follow the flow properly. 


3.2.1 .1 Tracer gas requirements 

The selection of a adequate luminescent gas is critical. Fundamental requirements 
that must be met are: 

- The luminescent lifetime must be long enough to trace the marked particles. Because 
a grid distortion of not more than 10% of the initial grid mesh size is requested, the 
required luminescent lifetime is a function of the expected velocities in the flow area 
of interest and of the initial grid mesh size. In order to obtain sufficient information 
for a jet formed by a 1 cm nozzle exit diameter, the maximum mesh size is set to 
about 5 mm. The jet core velocity to work with, is 17 m/s. Thus a luminescent life- 
time of at least 60 ps (twice the lifetime to detect the grid after 10% distortion) is re- 
quired. 

- The luminescent emission must be of sufficient magnitude to permit detection of the 
'distorted grid'. Luminescent intensity at a certain time after excitation, is a function 
of the quantum efficiency, the number of emitting molecules and the luminescent 
intensity decay over time. The latter, an important point for the tracer selection, has 
to be neglected because for most chemicals this function is not available. However, 
typically luminescence intensity decays exponentially over time. Thus, the relative 
luminescent intensity (normalized with the initial luminescent intensity) of a longer 
time emitting chemicals is expected to be higher compared to a shorter time emitting 
chemicals. If the laser intensity is high enough, the number of molecules emitting 
can be raised by seeding the flow with more luminescent molecules up to a limit that 
is given by the vapor pressure of the substance (at flow temperature). Thus a high 
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vapor pressure (commonly given at room temperature) of the luminescent substance 
is desired. 

Further criteria for the tracer selection which, although not critical, can ease 
experimental implementation: 

- to permit experiment implementation in air, luminescent emission should not be 
quenched by nitrogen, oxygen or other air components. 

- to ease illumination and detection problems, the absorption and emission wavelength 
should be close or in the visible region. 

- low cost 

- no toxicity 

The foregoing requirements are difficult to meet simultaneously. One drawback 
common to nearly all luminescent processes with long emission lifetime of chemicals in 
the gas phase, is the high quenching rate by oxygen. Also, substances with high vapor 
pressure tend to emit in the ultra violet, while those emitting in the visible tend to have 
vapor pressures considerably below 133.32 Pa (1 Torr). The most promising process is 
found in the phosphorescent emission of 2,3 butanedione (CH 3 COCOCH 3 ) also known as 

biacetyl. 

Because biacetyls phosphorescence is strongly quenched by oxygen, nitrogen is 
chosen as the carrier gas (see quenching rate constants in chapter 3.2. 1.2.2). 

The existing light source, a excimer laser which can emit light of 248, 308 and 
351 nm with different gas fillings, can excite biacetyl at the border of its phosphorescing 
first allowed absorption region (350-465 nm). This may reduce the achievable quantum 
yield and lifetime of biacetyls phosphorescence. If the direct excited phosphorescence of 
biacetyl can not be detected, biacetyls phosphorescence can also be sensitized efficiently 
by triplet sulfur dioxide, generated by excitation of sulfur dioxide within a band of 240- 
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320 nm. The sulfur dioxide sensitized phosphorescence of biacetyl by excitation of a 
nitrogen-biacetyl-sulfur dioxide mixture with a wavelength of 308nm is expected to have 
a higher quantum yield and lifetime than the direct excitation of biacetyl in a nitrogen-bi- 
acetyl mixture. 


3.2. 1.2. Biacetvl 

Biacetyl (2,3 butanedione, CH 3 COCOCH 3 or Ac 2 ), a watery and yellowish liquid 
with a strong, distinctive odor, is used as a flavor in the diary industry. It is nontoxic 
(average U.S. adult's daily consumption is about 10 mg), freezes at about -3 °C and boils 
at 89 °C. Its vapor is flammable with a flash point of 70 °C. Biacetyl is one of the few 
substances whose luminescence behavior is quite similar in all three phases, solid, liquid 
and vapor. Therefore, among other reasons, it has been extensively studied over the past 
fifty years. 

Biacetyl has at least three major advantages for use as a seed material for gases in 
laser marking: 

- It has a relatively high vapor pressure of 5.3 kPa (40 Torr) at 20°C and shows no 
obvious condensation below 2.66 kPa (20 Torr) at 20°C, so that a gas flow can be 
easily seeded with it 

- Its phosphorescence quantum yield O p (defined as emitted intensity/absorbed inten- 
sity) can be as high as 15 % [7,17,18,19] 

- Its phosphorescence lifetime can be reasonably long (practically up to 1.77 ms) 
[7,19] 

Furthermore, the relatively low triplet energy of biacetyl (Ef=57 kcal/mol) allows 

efficient triplet energy transfer from a variety of molecules of photochromical interest 
(like S0 2 ). 
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The major shortcoming of biacetyl is the severe affection of its phosphorescence 
quantum yield and lifetime by oxygen quenching (see figure 12). 


3 . 2. 1.2.1 Photochemistry of biacetvl vapor and biac e tvl vapor containing mixtures in a 
static cell 

The photochemistry of biacetyl vapor is one of the most thoroughly studied and 
best understood of all the carbonyl compounds. This favorable situation exists largely as 
a result of the many definitive studies by Noyes and his colleagues [20,21,22,23,24]. 
More recently, the direct excited as well as the sulfur dioxide sensitized phosphorescence 
of biacetyl vapor and biacetyl vapor containing mixtures has been studied by Horowitz, 
Calvert, Sidebottom, Kommandeur, Moss and coworkers basically in cell experiments 
[17, 1 8,25,26,27,28,29,30,3 1 ,32,33,34]. 

It shall be mentioned, that Horowitz and Calvert 1972 [17,26] found, that the 
observed effect of biacetyl pressure on the phosphorescent quantum yield in their cell 
experiments [26,35] resulted from biacetyl triplet decay at the cell wall, and concluded, 
that at least a major share of the earlier in cell experiments observed effect has the same 
origin. 


3.2, 1.2.2 Photochemistry of pure biacetv l vapor in a static rHI 

Biacetyl vapor displays a first allowed absorption band (350 - 465 nm) with a 
peak near 420 nm and a second allowed absorption band (220 - 320 nm) with a peak near 
275 nm. Whereas the direct excitation of biacetyl within its second allowed absorption 
region produces no detectable emission [18], excitation of biacetyl within its first 
allowed absorption band creates fluorescent emission, a singlet-singlet transition which 
occurs between states of the same multiplicity (i. e., states of the same electronic spin) 
and phosphorescent emission, a triplet-singlet transition which occurs between states of 
different multiplicity. 
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Figure 7: Absorption and emission spectra of biacetyl. The magnitude of absorption 
is measured by the extinction coefficient e. If and Ip are relative emission intensities 

at low pressure and room temperature for fluorescence and phosphorescence, respec- 
tively. Note the scale differences [201 


Fluorescence displays a emission range from 440 to 600 nm [20]. The fluorescent 
quantum yield of 0.25 % is essentially constant over a wide range of pressures (13.23 - 
5333 Pa), temperatures (25 - 119 °C) and exciting wavelengths (385 - 450 nm) 
[7,20,21,29]. The lifetime of the fluorescent emission is reported to be 10"^ to 10'^ s at 
25 °C for excitation from 365 to 435 nm and pressures from 40 to 5333 Pa [5,18,23,30]. 
Unfortunately no temperature dependent lifetime studies have been reported for the 
fluorescence of biacetyl vapor. The fluorescent emission of biacetyl vapor is not 
quenched by oxygen [30]. 

Biacetyls phosphorescent emission region extents from 490 nm to 600 nm with a 
peak near 510 nm. At room temperature, pumping near the absoiption peak of about 420 
nm and biacetyl pressures from 13 to 5332 Pa, biacetyl phosphorescence displays a life- 
time of 1.54 ± 0.23 ms. Its quantum yield of about 15 ± 3 % is essentially constant over a 
range of exciting wavelengths (365 - 440 nm) and pressures (1.5 - 533.32 Pa) 
[18,20,21,22]. At low wavelengths (< 380 nm) pressure must be high enough to prevent 
dissociation (> 267 Pa). 
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Unlike the fluorescence, the phosphorescence quantum yield and lifetime of bi- 
acetyl vapor is strongly dependent on temperature [17,18,19,23,27]. The phosphor- 
escence emission intensity as well as the lifetime decrease considerably with rising 
temperature above 50 °C. 



Figure 8 and 9: Logarithmic decay of biacetyl phosphorescence intensity and 
temperature dependence of biacetyl phosphorescence lifetime for 435 nm excitation. 
Data from Jian-Bang and coworkers [19] covers: total pressure 0.1 - 0.8 bar, biacetyl 
pressure 46 - 5507 Pa; [Ac 2 ]/[N 2 ] 0.001-0.30 


As a means of understanding the biacetyl luminescence process, it is useful to 
examine the excitation scheme (figure 10). 

The luminescent processes are as followed. Molecules can be pumped by absorp- 
tion of energy from the ground singlet state S 0 (’A g ) to the singlet state S ( 1 A ll ; zeroth 

level ca. 64 kcal/mol above the ground state [29]). Most of the molecules in the S state 
radiationlessly transit to the triplet T ( 3 A„; zeroth level ca. 57 kcal/mol above the ground 

state [23,26]) (quantum yield = 1 [7,22,29]) under the perturbation of the spin-orbit cou- 
pling (S — » T intersystem crossing, an intramolecular process and independent of pres- 
sure [7,22,29]), the other return to the S 0 state via either the S — ► S 0 internal conversion 

or the spontaneous transition emitting fluorescence. Some of the molecules in the T state 
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Figure 10: Biacetyl excitation scheme. f7,36] 


retum to the ground state S 0 via T — > S 0 intersystem crossing or the phosphorescence 
transition (quantum mechanically forbidden with a long radiative lifetime of 10 ms) with 
a actual observed lifetime of 1.54 ± 0.23 ms at 25°C [18,29], The other either transit to S 
via the T — » S intersystem crossing or dissociate into 2(CH 3 CO) (about 70 kcal/mol for 

the ground state molecule necessary to dissociate [26,29]). It is also possible that the 
biacetyl molecules in the T state are quenched by colliding with molecules of other 
species, but except a few species including the 0 2 , the quenching rate constant is quite 

small [17,18,19,29]. 

The actually observed lifetime of the phosphorescence transition is shorter than 
the radiative lifetime due to two important processes: 

A annihilation reaction between two molecules in the triplet state (quenching rate 
constant: k=4.2 ± 1.7 X 10 u cm 3 /(mole*s) [31]). The significance of this path in- 
creases with concentration of triplet molecules and, hence, also with laser power. 
It should be emphasized, however, that although the time constant is shorter at 
higher laser intensities, the absolute phosphorescence yield is still enhanced, if the 
laser intensity used does not reach the saturation value of triplet biacetyl forma- 
tion. Thus, high laser intensities might be desirable [7,18]. 





-25- 


Collisional de-excitation (quenching) also shortens the observed phosphorescence 
lifetime. Self -quenching, quenching by sulfur dioxide and quenching by nitrogen 

are far less important than quenching by oxygen (quenching rate constants in 
cm 3 /(mole*s): k^s 4 X 10* [18]; k S02 = 4 X 10* [18]; k N2 = 10* [7]; k^-5 X 

10 11 [17,18]). The latter is a strong quencher essentially precluding use of 
biacetyl in systems containing air. On the other hand, the low nitrogen quenching 
rate allows operation over a wide range of nitrogen background pressure, a fact 
that was verified in a static cell at total pressures between 1.3 kPa (10 Torr) and 
100 kPa (1 atm). [7,36] 


3.2. 1.2.3 The S0 7 -sensitized Phosphorescence of bi a cetvl vapor in a static cell 

Sulfur Dioxide has three main regions of absoiption in the near ultraviolet: A first 
very weak forbidden absoiption band from 340 - 390 nm, a stronger, first allowed 
absorption region from 240 - 320 nm and a still stronger absoiption in the 190 - 220 nm 
region [28,37,38,39]. 

Static cell experiments by Horowitz, Calvert and coworkers [24,25,26,34] have 
demonstrated that the excited triplet state of sulfur dioxide ( 3 S0 2 ) is the primary reactive 

entity formed in the photolysis of sulfur dioxide irradiated within the second allowed 
absorption band (240 - 320 nm). Die main route by which triplet sulfur dioxide 
molecules are formed is a second-order intersystem crossing reaction involving the sin- 
glet sulfur dioxide molecule (’S0 2 ) and some collision partner. 

Within the first allowed absorption band of S0 2 (240 - 320 nm) irradiated mix- 
tures of sulfur dioxide and biacetyl demonstrated, that the incoming light was only ab- 
sorbed by S0 2 , and that the phosphorescence in biacetyl can be sensitized by triplet sul- 
fur dioxide molecules. 


26- 


Even small amounts of biacetyl 
added to within its first absorption band 
irradiated sulfur dioxide show a decrease 
in the phosphorescence of sulfur dioxide 
(380 - 470 nm) and a sensitized quantum 
yield of biacetyl emission Rao 

and coworkers [27] showed, that relative 
small biacetyl pressures (5 Pa) quenched 
the phosphorescence (lifetime = 0.4 ms 
[28,37]) of sulfur dioxide (84 - 293 Pa) 
completely and did not lower detectably 
the emitted fluorescence (lifetime = 

0.02ms [27]) of sulfur dioxide. Note in 
figure 11, that in the mixture the 
phosphorescence of S0 2 is quenched 

considerably and the sensitized phos- 
phorescence emission of biacetyl appears (quenching rate constant for 3 S0 2 quenching by 
biacetyl: 1.42 X 10 -11 l/(mol*s)). 

For [SOJ / [AcJ concentrations from 60 to 600 and S0 2 - pressures of 85 - 550 
Pa, Horrowitz and Calvert [25] found in static cell experiments, that the reciprocal of 
the quantum yield of the sensitized phosphorescence in biacetyl is, respectively, a linear 
function of the reciprocal of the total pressure (for [SOJ / [Ac 2 ] = const.) and of the 
concentration ratio (for S0 2 - pressure = const.). 

Whereas the dependence on the [SOJ/fAcJ ratio is anticipated theoreti- 
cally, the dependence on the total pressure for fixed concentration ratios is 
theoretically unexpected and is probably largely the result of biacetyl triplet diffusion 
with deactivation at the cell wall (diffusional effects at the cellwall decrease with pres- 
sure) [17,26]. Unfortunately new data from experiments in larger cells has not been re- 
ported. 
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Figure 11: Typical low-resolution spec- 
trofluorimeter trace of the fluorescence 
and phosphorescence emission from pure 
sulfur dioxide and from a mixture of S0 2 

and biacetyl [271. 
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The lifetime of triplet bi- 
acetyl phosphorescence in 
biacetyl-S0 2 and biacetyl-oxygen 

mixtures has been investigated by 
Sidebottom and coworkers [18], 
by exciting the mixtures within 
the second allowed absorption 
band of biacetyl. They showed, 
that the lifetime of biacetyl triplet 
in biacetyl-S0 2 mixtures is rela- 
tively insensitive to change in 
added S0 2 pressure for biacetyl 

pressures of 1.6 Pa and 386 Pa, 
whereas biacetyl triplet is 
quenched very effectively even at small oxygen pressure. 

No lifetime measurements of biacetyl phosphorescence in within the first al- 
lowed absorption band of sulfur dioxide excited biacetyl-S0 2 mixtures have been pub- 
lished. 

3.2. 1.2.4 The ’excess' S0 2 -sensitized biace tvl phosphorescence at high added gas pres- 
sures in a static cell 

In static cell experiments Horrowitz and coworkers [25,38] determined the quan- 
tum yields of sensitized biacetyl phosphorescence emission for S0 2 -biacetyl-N 2 mixtures 
irradiated within the first absorption band of S0 2 . Experiments at constant concentration 
ratios ([N 2 ] : [AcJ : [SOJ = 255 : 1 : 1) with N 2 pressures from 0.08 bar to 0.7 bar at 
25°C were performed. 

The data, that can not be explained by the investigated mechanisms at low 
pressure (<0.01 3bar) [26,35], reveals the continues to rise as the total pressure of 

the gas mixture is increased fexcess’ biacetyl phosphorescence). Since there is no 
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Figure 12: Stem-Volmer plots of the reciprocal 
of the biacetyl triplet lifetime vs. the pressure of 
added oxygen and sulfur dioxide gas. [181 
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significant increase in quantum yield of 
phosphorescence on increasing the pressure 
of added N 2 gas from 0.01 to 1 bar in bi- 
acetyl-nitrogen mixture, a S0 2 species 
seems to be responsible for the 'excess' bi- 
acetyl phosphorescence. 

The results reported in Horrowitz and 
coworkers work [25,38] offer support for the 
existence of a undefined species X. In irradi- 
ated S0 2 systems at high pressure it is sug- 
gested, that X is not a triplet species which 
can transfer energy direcdy to biacetyl, but 
may form 3 S0 2 rather efficiently on collision 

with added gases. 

Unfortunately no lifetime studies have been reported for the 'excess' sulfur diox- 
ide sensitized phosphorescence of biacetyl at high added gas pressures. Since the direct 
excited phosphorescence lifetime of biacetyl is, respectively, not influenced by high 
added N 2 pressure and added S0 2 pressure, the lifetime of the S0 2 -sensitized 
phosphorescence of biacetyl is be expected to be independent on added N 2 or S0 2 (above 
a minimum S0 2 pressure) pressure. 
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Figure 13: Plot of the reciprocal of 
the quantum yield vs. the reciprocal 

of the total pressure for data from a 
N 2 -Ac 2 -SQ 2 mixture photolysis [25] 
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3.3 Performed experiments 

LIPA experiments and the lifetime experiments are performed with nitrogen-bi- 
acetyl gas mixtures ([N 2 ]:[Ac 2 ]=200:3) and the set-up that is shown in figure4. The 

experiments with a still camera, that replaces the ICCD-camera assembly in figure 4, are 
performed with varying concentration ratios of nitrogen-biacetyl-sulfur dioxide gas mix- 
tures and different flow generating set-ups (see figure 14). The core jet velocity for all 
experiments is, respectively, 17 m/s. A co-axial annular jet is provided if the main jet is 
issued into a tube. Relating to a preliminary experimental set-up and space, the annual jet 
is supplied further away from the initial jet as desired. The gas mixture component 
concentration in the side jet is, respectively, the same as in the main jet. 

The concentration ratio for biacetyl in the gas mixture is measured by the de- 
crease of liquid biacetyl in the mixing chamber during a certain time interval. The sulfur 
dioxide flow ratio is measured by a flowmeter in the sulfur dioxide support line. 


3.3.1 Still camera experiments 

In order to optimize the detectable biacetyl emission, experiments using a still 
35mm camera instead of the ICCD camera (see figure 4) are performed. A uniform grid 
pattern brightness over a wide area in- and outside the potential jet core is desired. The 
camera is equipped with a 58mm (f-stop: 1:1.2) and a close-up lens and is loaded with 
3200 ASA (Kodak T-MAX P3200) film material. The camera-shutter and the laser are 
operated manually. To capture all the emitted light, the camera-shutter is opened before 
the laser is fired and closed after the emission disappears. Pictures of the whole grid, as 
well as pictures of only one line (blocking of the other grid-lines, hence same energy in 
this one line as in one line of the hole grid) are taken. To keep the experimental set-up 
simple, we tried to find a set-up that allows the experiment to be performed in an open 
system, although biacetyls phosphorescent emission is strongly quenched by oxygen. 
Five different set-ups that are shown in figure 14 were examined. 
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Figure 14: Schematic diagrams of the five examined flow field set-ups. See appendix 

F for more details. 







Using a Pitot-tube, the main jet velocity is measured to be .respectively, 17 m/s in 
the potential core. Hence, performing experiments with one grid-line, considering the 
expected flow-field and the lifetime (1.8 ms) of the sensitized biacetyl emission, pictures 
taken of the emission are sumilar to that shown in figure 15. Whereas the excited 
particles in the central core of the jet are supposed to move with decreasing emission 
intensity (see chapter 3.2. 1.2.2) over a distance, that is given by the detectable lifetime of 
the emission and by the jet velocity, the distance excited molecules in the mixing area are 
expected to move decreases with the distance from the central core because the mean 
velocity decreases. 

Nitrogen gas as well as nitrogen-sulfur dioxide, nitrogen-biacetyl and nitrogen-bi- 
acetyl-sulfur dioxide gas mixtures of the initial jet (same gas mixture in side jet) have 
been examined to determine the performance of the sensitized biacetyl phosphorescence. 
Different sulfur dioxide concentrations at constant biacetyl concentrations in the initial 
jet have been examined in order to optimize the detectability of the emission. 
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The different flow generat- 
ing set-ups (see figure 14) have 
been employed to reduce the 
oxygen concentration in the 
measuring area. Lower oxygen 
concentration, hence better 
detectability because of less 
quenching of the phosphorescent 
emission, is expected for issuing the 
jet into a tube (set-up: 2). A longer 
tube-length (set-up: 3) and a 

pressure gradient providing 
honeycomb in the tube outlet (set- 
up: 4) is expected to decrease 
oxygen concentration further. The 
pockets (set-up: 5; see appendix F) 
are added to decrease the energy 
loss along the incoming laser lines 
before reaching the measuring area 
(see chapter 4.1). 

The taken pictures are presented in appendix B. 
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Figure 15: Expected pictures (black area) taken 
by a still camera schematically. 


13.2 ICCD camera ex periments 
3.3.2. 1 LIPA experiments 


For the data sets presented here, an overall measurement area of approximately 
the distance from the jet centerline of ±3 y/r, and the distance from the jet nozzle in 
stream wise direction of 0 to 8.5 x/r was used. 18 incoming laserlines create 47 
intersections and 30 grid boxes in the measuring area. The mesh size of the generated 
grid boxes spreads from about 4.2 to 5.2 mm 
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Figure 16: Four raw LIPA data grids of a 17m/s nitrogen-biacetyl jet, taken with 
delay times of 0, 30, 60 and 90 us after the laser pulse. 


The chosen measuring area provides a data resolution of 0.0735 mm/pixel on the 
detector. The data frames are taken with a gating time of 0.095 ps, providing a maximum 
line movement during gating of a distance smaller than 1.7 pm, hence a maximum 
movement of 2.5% of the associated pixel distance on the detector. Whereas reference 
grids are taken when the laser fires, distorted grids are take at 30, 60 and 90 ps after the 
laser pulse (four raw data frames are shown in figure 16). 


The excimer laser is set to provide 220mJ during a pulse length of 20ns with a 
repetition rate of 3 pulses per second (internal laser trigger). Each laser line carries 
approximately 6.6mJ. The grid plane is located above the centerline of the jet exit nozzle 
and has a thickness of the same order as the width of each beam (about 0.5mm). 

Because the emission of a by 308nm excited nitrogen-biacetyl mixture of a 
molecular concentration [N 2 ]:[Ac2]=200:3 ( mass concentration of biacetyl - 5 %) can be 

easily detected with uniform brightness over the measurement area, the data grids are 
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taken without added sulfur dioxide gas. Experiments made with added sulfur dioxide gas 
show a ununiformity of emission over the measured area. The brightness of the lines in 
the potential core area of the jet increased more than in the surrounding area. Adjusting 
the camera intensification for the bright lines involve less contrast in the darker areas, 
hence less detectability in the darker area. To take maximum advantage of the dynamic 
detector range, a uniform line brightness in the hole picture area is required. 

Pockets, as they are employed in the still camera experiments to decrease the 
energy loss along the laserline, are not employed in the LIPA experiments. The 
disturbance of the flowfield is considered to be not negligible, hence no pockets are used. 


3.3.2.2 Lifetime experiments 

In order to measure the detectable lifetime of the biacetyl emission in a set-up, 
that is the same as in the LIPA experiments, the ICCD camera is set to maximum 
intensification. The same gain duration as in the LIPA experiments (0.095jjs) is 
employed. To generate a higher contrast in the picture area, the side-jet is not seeded 
with biacetyl. In order to gain a even higher signal-noise performance of the expected 
weak signals, 5 pixels along both axes are software binned. 

Binning in software is chosen because the hardware binning uses shift registers 
and binning capacitors that are only approximately twice the size of regular pixels. In 
order not to loose the true signal value, hardware binning is therefore limited to 2 pixels. 
Implementing software binning, the CCD is read out using combined pixels sufficiently 
small to avoid saturating the shift registers or binning capacitors. These stripes are then 
added in the software producing a resultant combined pixel that represents more photons 
than is possible only using hardware binning. 

The signal noise performance is generally improved by binning because the extra 
noise collected with each additional stripe read is negligible compared with improvement 
in the photon signal to noise ratio due to the overall large number of photons collected. 
Although binning reduces the image resolution, wherefore it is not been implemented for 
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the LIP A experiments, it is applied for the lifetime experiments where high 
intensification and signal to noise ration is more important than high resolution. 

Examples of pictures taken are shown in appendix D. To prevent the CCD array 
from overexposure, pictures for time delays smaller than At=0.5ms are not taken. Smaller 
delays are associated with brighter phosphorescence. 
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4. Results 

4.1 Still camera experiments 

Unlike experiments using nitrogen-biacetyl-sulfiir dioxide mixtures 
([N 2 ]:[Ac 2 ]:[SO 2 ]=200:3:2 to 200:3:5), pure nitrogen gas, nitrogen-sulfur dioxide 
([NJ : [SC>2]=200: 1 to 200:5) and nitrogen-biacetyl ([N2]:[Ac2]=200:5) gas mixtures of 
the initial jet do not emit detectable emission, whichever flow generating set-up (see 
appendix B) is used. The maximum length of detectable emission in stream wise direction 
is approximately ^.=25mm. Considering a core velocity of 17m/s, a 'lifetime' of about 

1.5ms can be detected with the applied equipment (see chapter 3.3.1). Because of the 
detected 'lifetime' and the absence of detectable emission in nitrogen, nitrogen-biacetyl 
and nitrogen-sulfur dioxide gas mixtures, the detected emission in experiments using 
nitrogen-biacetyl-sulfur dioxide gas mixtures is assumed to origin by sulfur dioxide 
sensitized phosphorescence of biacetyl. 

Performed experiments with set-up: 1 show that the emission in the central core 
of the jet can be easily detected, whereas there is no detectable emission in and outside 
the mixing area of the jet even with high present sulfur dioxide concentrations (up to: 
[SC>2]=1.3 X [AcJ) in the initial jet. The detectable emission increases by rising the 
sulfur dioxide concentration from [N 2 ]:[Ac 2 ]:[S0 2 ]=200:3:l to 200:3:5. Quenching of 
the biacetyl phosphorescence by oxygen is assumed causing this lack of detectable 
emission in- and outside the mixing area. 

Issuing the jet into a short tube (set-up 2; tubelength: 180mm) increases the 
detectable emission in the central core region of the jet reasonably. The same tendency of 
emission intensity with increasing sulfur dioxide concentration as in set-up 1 can be 
observed. However, outside of the jet area no emission can be detected. Quenching of the 
biacetyl phosphorescence by oxygen because of the relative long distance between the 
side jet and the main jet is assumed to be responsible for this absence of detectable 
emission outside the jet. 
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By issuing the jet into a longer tube (set-up 3; tubelength: 740mm), the intensity 
of the detectable emission in the central core of the jet reduces to a level less than for set- 
up 2, however, emission outside the potential core can be easily detected. The absorption 
along the incoming laser line, outside the potential core, is made responsible for the 
reduced emission level in the central jet area. Because in set-up 3 emission even outside 
the jet can be detected, the absoiption along the laser line in- and outside the shown area 
may be reasonably higher than in set-up 2, hence less energy in the incoming laser line 
and less detectable emission. 

A maximum of emission intensity between [NJrfAcJrfSOJ ratios of 200:3:1 and 
200:3:3 is detected for experiments with set-up 3. The picture area in appendix B shows 
about 30% of the total laser line inside the tube. Assuming homogeneous gas 
concentrations along the incoming laser line and no saturation effects, the decay of the 
laser intensity Ij along the laser line [41] coordinate x, may be approximated by 

Il~Il 0 *exp(-ax,). Assuming the absoiption coefficient to be linearly dependent on the 
sulfur dioxide concentration in the illuminated gas volume a=C 2 *[S0 2 ] and a linear 

dependence of the emitted light intensity on the sulfur dioxide concentration and the 
illuminating light intensity I e ss Ci*Ij*[S 02 ], the emitted light intensity I e may be 
roughly approximated by the relation I e =C 1 *Ii 0 *[SO 2 ]*exp(-C 2 *[SO 2 ]*x 1 ). This 

relation shows the existence of a maximum emission intensity at a constant position on 
the laser line as a function of sulfur dioxide concentration in the gas mixture. However, a 
maximum emission intensity may not occur in the experiments if the optimum sulfur 
dioxide concentration for a x-coordinate in the picture area is not exceeded. Experiments 

with set up 1 and 2 do not show a existence of a maximum emission intensity because 
the origin of the coordinate x, is at the edge of the jet, hence lower x r values. 

Nevertheless this may explain the existence of a maximum emission intensity as it is 
observed in experiments with set-up 3. 

A further decrease of detectable emission of excited nitrogen-biacetyl-sulfur 
dioxide mixtures occurs by plugging the tube outlet of set-up 3 with a honeycomb. An 
overall decrease of detectable emission in comparison to set-up 3 is observed. As well as 
in the experiments using set-up 3, a maximum of the emission intensity can be detected 
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within the applied concentration ratios in set-up 4. The same explanation as in the 
foregoing paragraph shall be given. 

In set-up 5 'pockets'(see appendix F) arc implemented to rise the laser line energy 
in the picture area, by reducing the length on which absorption of the laser line energy 
can occur before entering the picture area by about 50%. Obviously the detectable 
emission intensity is much higher than in experiments with set-up 4, which supports the 
above given explanations. 


4.2 LIPA e xperiments 

From initially 40 recorded data frames, 24 data frames of distorted grids with 47 
intersection points each are reduced on the image processor. The data frames of the 
distorted grids are taken 60ps after the reference grid which is, respectively, taken at the 
time of the laser pulse. The gaiting time for the camera intensifier is adjusted to 0.095ps. 
The fluid dynamic properties which are the mean velocity components, the Reynolds 
stress, the turbulent intensity and the vorticity are calculated. In appendix C the results 
are shown in form of tables, vectorplots and contourplots. 

Because the data reduction on the image processor is performed by hand (finding 
the grid intersection points with a software cross on a 1024x1024 pixel screen) a random 
error in locating the intersection points occurs. This uncertainty in location the proper 
intersection points is assumed to generate the major error in the experimental data and is 
therefore particularly examined (see also chapter 3.1.3). In order to evaluate the reading 
error and hence the experimental error, several readings of different intersection points in 
a variety of different data frames are performed. The maximum deviation from the 
average location of the readings is determined to 3 pixel of the initial data resolution for 
the distorted grids and to 1 pixel of he initial data resolution for the reference grid. 

The difference in the position of the distorted and undistorted intersection points 
occurs because of the fluid motion during luminescent emission decay with time and 
because of diffusion. Typically phosphorescent emission decays exponential with time 
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(see chapter 3.2. 1.2.2). Applying the same dynamic range (14-bit detector, but just 8-bit 
for data reduction) for different maximum light emission intensities at unchanged 
background light intensity will reduce the contrast, the background noise or both. 
Because a small contrast as well as a high background noise render reading more 
difficult, hence, less reproducible and accurate. The broadening of the visible lines, 
which happens as a result of diffusion (mainly turbulent diffusion) that on average carries 
emitting fluid volume out of the excited fluid volume does not necessarily influence the 
reading accuracy of the intersection points (the reading of the center point is still accurate 

as long as the border of the area can be well determined), but the associated decrease of 
contrast reduces accuracy. 

The measured mean velocity distribution is shown in appendix C. Velocities of 
the order of the expected core velocity are measured at measuring points in the area 
around the jet axis (ly/rl<1.2). Slow, negative streamwise velocity is determined for 
measuring points in the interval 2.2<ly/ri<2.9, whereas slow streamwise velocity is 
determined for the other measuring points. The measured y-component of the mean 
velocity is only larger than lm/s at some measuring points at ly/rl>2.2. 

Considering a circular jet of uniform velocity coming out of a nozzle into a large 
stagnant mass of the same fluid, the generated flow field is commonly sectioned into a 
potential core, an annual shear layer and the ambient fluid. Due to the velocity 
discontinuity at the plane of the nozzle, shear stresses are set up and a shear layer 
originates. In most of the practical cases, this shear layer becomes turbulent very close to 
the nozzle itself. On the inside it invades into the uniform velocity jet and penetrates into 
the ambient fluid on the outer side. 

In order to compare the measured velocity field to other experimental results in 
the literature, the measuring area is divided into the three common areas by the empirical 
equations [42]: 


-39- 


|^| = 0.95 - 0 . 097 «^ ( 4 . 1 ) 

|fJ = l- 07 + 0.158 * 7 (“. 2) 

with: y, - inner edge of mixing region 
y 2 - outer edge of mixing region 

Whereas the core velocity (inside the inner edge of the mixing region) is expected 
to be constant and equal to U c , the velocity outside the mixing region is expected to be 

equal to zero. Inside the mixing region, the velocity distribution is predicted with the 
cosine function of Squire and Trouncer [42]: 



(y 2 -yV(y 2 /yi) 


Figure 17: Comparison of the measured velocity to the empirical equation of Squire 
and Trouncer [42]. ▼: here performed measurements with rms-eiTors; — : Squire and 
Trouncer (equation 4.1 to 4.3). 


Denoting the average of the streamwise velocities in the estimated central core as 
the core velocity U c leads to the distribution in figure 17 and to the deviations (U 

TT \ U • . V Xcst * nit!e 

J shown in table Cl. Inside the central core region velocity deviations up to 1.28m/s 

can be observed, whereas maximum velocity deviations of 4.79m/s and 5.23m/s occur in 
the mixing and ambient region. 
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Agreement with the predicted velocity distribution in the potential core region 
and in the inner half of the mixing region is better than in the area further away from the 
jet-axis. The observed deviations from the empirical velocity distribution as well as the 
measured negative stream wise velocities in the ambient region imply a significant 
influence of the surrounding tube and the secondary flow on the velocity field. 
Particularly the measured negative stream wise velocities hint to a recirculation zone 
typically observed for jet mixing in a duct [10]. 

The general result of high velocity jet mixing in a duct is increasing static 
pressure along the flow direction. The produced pressure gradient can effect a massive 
rearrangement of the flow. Selecting a variable area duct to keep the static pressure 
constant can prevent recirculation. Here, however, an axial pressure gradient is required 
to obtain a low oxygen concentration in the measuring area. Special conditions have been 
determined in the literature where the flow remains 'similar’ even in the presence of an 
axial pressure gradient [10,43,44]. 

The velocity, turbulence intensity and Reynolds stress data obtained have been 
compared to measurements by Hussain and Clark [45], Crow and Champagne [46] and 
Sami, Carmody and Rouse [47]. The Mach numbers of the free jets in the reference- 
experiments is sufficiently low that the flow can be considered incompressible. The 
Reynolds numbers of these experiments are higher than in the present measurements 
(^ e d — 360000, 106000 and 220000 compared to 12000), but lower Reynolds number 

data for the developing region of a free jet has not been found in the literature. In order 
to compare the measured fluid dynamic quantities to the literature-data, the measuring 
area is sectioned into 7 intervals in flow direction. Plots of the quantities in this 7 
intervals are shown in appendix C (figure Cl 1). 

A peak of the Reynolds stress and turbulent intensity distribution at ly/rl=l as well 
as a minimum at ly/rl=0 is predicted [9,10,42,45,46,47]. Neither of these predictions can 
be observed in the actual data. Compared to the calculated Reynolds stress, the calculated 
turbulent intensity is rather equally distributed over the measuring area. Both do not 
show a significant tendency or symmetry. 
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As well as the Reynolds stress and the turbulent intensity, the mean vorticity 
distribution is expected to have a maximum at ly/rtel and a minimum at ly/H=0. n, e 
latter can be easily observed in figure C9 and CIO. For the 6 measuring points close to 
e jet symmetry-axis the measured vorticity expands from -241 to 340 1/s. The vorticity 
seems to increase with the distance from the jet axis to a maximum between l<ly/ti<2. 
The calculated vorticity-data shows qualitative as well as quantitative symmetry to the jet 


4,3 Lifetime exnerimei^ 


Because the slow co-axial jet that formed an annulus jet has not been seeded with 
■acetyl, only the potential core and the mixing region of the 1cm central jet can be seen 
as the bright area in the pictures shown in appendix D. 


The jet area can be easily observed with decreasing light intensity up to a delay 
rime of 2.5ms. For longer time delays the jet contour disappears slowly and the light- 
intensmes of the background and the jet become similar. For the implemented flow 

condmons and detection equipment, a reasonable limit of dme delays applicable for 
LIPA is found to be 2.0 to 2.5 ms. 


The illumination of the jet is performed with the same laser grid pattern as in the 
IPA experiment. The impressed grid pattern can not be observed. Two factors can be 
thought of being responsible for the disappearance of the grid pattern: 

- The resolution of the pictures made here is lower than in the LIPA-experiments. 
5x5 pixels of the resolution implemented in the LIPA experiment are binned (see 
chapter 3.3.2.2). With this resolution a initial grid-box has a sidelength of 
approximately 15 binned pixel blocks. This may cause likely smoothing of the 
image, but can not be made individually responsible for the lack of grid pattern. 

Diffusion causes exchange of emitting and non emitting molecules between the 
excited and not excited flow volumes. By the exchange of the molecules the 
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impressed lines widen. Because the diffusion is a function of time, the width of 
the emitting lines increases with the time delay At The diffusion length l d of the 

excited molecules during the delay time At=60ps can be calculated with the 
solution for the boundary layer of the suddenly accelerated plane wall ([48], page 
83): 


Id = V v* At = 0.03 mm (4.4) 

Compared with the full image resolution of 0.0735 mm/pixel (with 
binning only 0.3675 mm/pixelblock) the diffusion should be negligible as well as 
the maximum line-movement of 0.0016 mm during the 95 ns gating time. The 
magnitude of the calculated diffusion length does not support considerable 
influence of the diffusion on the disappearance of the grid pattern. 
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5. Discussion 

If illumination, hence grid generation, shall be performed with a wavelength of 
308nm, the use of sulfur dioxide to sensitize biacetyls phosphorescence might be 
desirable for the implementation of the LIPA-technique in a gas flow. As in the static cell 
experiments reported in the literature and in the presently employed experimental 
conditions, the sulfur dioxide-concentration in a 308nm irradiated sulfur dioxide- 
nitrogen-biacetyl mixture can effect the emitted light intensity significantly. 

For the present investigation, the stimulation process and concentration ratios of 
the mixture components, the presence of sulfur dioxide increased irradiated light 
intensity considerably. However sulfur dioxide can increase the efficiency of the 
stimulation, the use of the sensitizer sulfur dioxide for the LIPA-technique requires 
adjustments that may not always be matched. The observed decrease of emission 
intensity along a laserline even for small sulfur dioxide concentrations, makes a 
adjustment of mainly two parameters necessary. The length, laser beams pass through the 
test atmosphere, and the sulfur dioxide concentration in the gas mixture have to be 
adapted. 

In principle, absorption along the laser line always takes place if excitation 
occurs. Direct excitation of biacetyl closer to its absorption peak by a dye-laser may 
cause a similar effect as remarked. Higher absorption and hence higher decrease of 
energy may cause observable decrease of emission intensity along a incoming laser beam 
as a function of biacetyl pressure. Because implementation of the LIPA-technique desires 
homogeneous emission intensity of the laser lines in the measuring area, the different 
experimental parameters are to be optimized. In general a high quantum yield, hence a 
high conversion rate of absorbed to emitted intensity, is desired in order to generate 
bright and homogeneous grid-lines. 

The present performed LIPA-experiments are performed with a nitrogen-biacetyl 
gas mixture. Sulfur dioxide is not employed because inhomogeneous emission intensities 
in the measuring area are observed by adding the sensitizer of biacetyls phosphorescence 
to the gas mixture. Besides the described absorption along the grid-lines, mainly 
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mhomogeneous gas mixture over the measuring area is made responsible for this effect. 
Considering, respectively, the similar molecular weights of biacetyl and sulfur dioxide 
the concentration distribution of this two gases is assumed to be similar over the! 
measuring area. The concentration of these two gases is expected to be the highest in the 
central core, where the oxygen concentration is supposed to be the lowest. Support for 
this assumption is provided by the high emission intensities in the central core that could 
have been observed within the still camera experiments. Taking the limited dynamic 
range of the detector and the high differences of emission intensities over the measuring 
area in account, can explain the inhomogeneity of the detected emission intensity A 
sulfur dioxide concentration that provides reasonable signal-noise improvement without 
generating the mentioned disadvantages has not been not found for the employed set-up. 

fso^rMi^m 0 / T Cra experiments U is to be in the range of 

fSO 2 ]<[N 2 ]/200 for [N2]:[Ac2]=200:3. However, a sulfur dioxide-concentration ratio in 

this range can not be properly adjusted in the here employed set-up. 


The presented error analysis (chapter 3.1.3) permits evaluating the accuracy of 
the calculated fluid dynamic properties. The error analysis predicts a tendency of 
decreasing accuracy in the obtained velocity, vorticity, turbulent intensity and Reynolds 
stress data which is supported by the experimental results. The plots in appendix C show 
a qualitative as well as quantitative high rate of symmetry for the velocity and vorticity 
field. In contrast no general tendency of the turbulence intensity and Reynolds stress data 

could be observed. Furthermore, a random Reynolds stress distribution with a peak that 
is not explainable is observed. 


In order to reevaluate the accuracy of the obtained data, the reading error is 
approximated by the observed velocity deflection in the central core (AU x =1.28m/s; see 

chapter 4.2). Neglecting errors of other origin than reading errors 8s 1 Ld 8s", htnct 
calculating reading errors by the obtained velocity deviation in the predicted central core, 
leads to a k posteriori evaluation of the obtained data. This appraisal guides to a 
deterioration of the & priori estimation of the data precision by a factor of about 5 
(AU x -1.28m/s k priori compared to the k posteriori calculated confidence region of 
0.25m/s). 
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A comprehensive explanation for this discrepancy can not be given. One 
important factor may be the duct effect on the flow field. Unsymmetrical, oscillating 
recirculation zones have been observed in a duct under some flow conditions by Curtet 
[10]. If a similar flow-field has been generated in this work, the observed velocity 
defection m the potential core may be generated by a oscillating potential core. 

Some suggestions for future improvement of the performance of the LIPA- 
techmque in general, and measurements in a jet in particular shall be made here: 

To resolve the velocity more uniformly over the whole measuring area, different 
time delays for areas of different velocity have to be implemented. If supported by 
the sampling rate of the equipment, this can be achieved by taking pictures of 
distorted grids at different time delays during one illumination process. If the 
equipment does not provide a fast enough sampling rate the areas of different 
velocity should be measured at separate excitation processes. The time delay At 
should be adapted to the velocities in the flow regions. 

In order to improve local resolution, smaller mesh sizes or a larger orifice diameter 
might be desirable. The creation of smaller grid sizes in the order of 1mm or even 
smaller, may require implementation of a different beam dividing technique. 
Easing of handling and adjusting as well as more distinguished grid-lines should 
be the main goals. The use of diffraction gratings or fiber optics may be worth to 
be further examined. A larger orifice diameter might be required to resolve the 

expected peaks in the fluctuation and vorticity distribution as they are shown in 
appendix C. 


Improvements in the flow generation have to be achieved. Since oxygen has to be 
kept out of the measuring area, an open air system will always require a pressure 
gradient. According to literature free jet characteristics as well as a pressure 
gradient can be accomplished in a duct. If the relative simple set-up of a open 
system shall be retained, changes on the duct have to be made to realize free jet 
characteristics. 
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- Implementation of a closed system for the flow generation. Issuing the jet into, or 
generating the jet inside a large, with nitrogen and biacetyl vapor filled, box can 
provide free jet characteristics. Generating the flow-field of interest inside a closed 
x would make a more precise adjustment of the gas mixing ratios as well as a 
more homogeneous mixing possible. Generating of the flow-field could be 
established by replacing the bottom plate of the tank used in the present 
experiments by a fan, that forces the flow through the nozzle. Placing this set-up 
mto a closed loop-system filled with the desired gas mixture is a suggestion that 

would also allow forestalling quenching of biacetyls phosphorescence by oxygen 
without the need of a duct. 

A more homogeneous gas mixture of nitrogen and biacetyl may be desirable. 

eahzing seeding with a carburetor or a fuel injection like functioning system can 
improve homogeneity as well as provide more accurate concentration ratio data. 

If a better evaluation of the grid-intersections is attained, a higher resolution of the 
image capturing device can improve accuracy of the technique. Film material in 
combination with image intensification or higher resolved CCD arrays may be 

helpful Considering automatization of the technique, a digital data acquisition 
rather than a analog technique should be preferred. 

Achieving more accuracy in the calculated fluid dynamic quantities requires 
uation of a much larger number of frames than it is done in this work, 
erefore automated data reduction is necessary. Computer programs that are 
currently being worked on may offer this in the near future. 
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6. Conclusions 

The experiments are conducted in an completely gaseous environment with 
properties very near to those of air. 

A jet in a recirculating co-axial flowing environment is studied. 

- Recirculation outside the generated jet can be observed. Re-entrainment of the 
central jet is likely. In order to provide complete free jet characteristic the set-up 
has to be changed. Changes on the duct according to literature or the use of a 
closed system are suggested. 

- On average, the velocity vector field picture for the predicted jet area of the 
central core and the mixing region looks like expected. High velocities in the 
central core and decreasing velocities in the mixing zone are measured. 

The tendency of the mean vorticity distribution perpendicular to the jet axis looks 
like expected. Small vorticity close to the jet axis, increasing positive vorticity in 
Positive y-direction and increasing negative vorticity in the negative y-direction. 
Unfortunately only one vorticity measuring point is far enough away from the jet 
axis to show the increase of vorticity for high y/r (>2). 

- The accuracy of the performed data reduction does not provide reasonable results 
for Reynolds stress and turbulence intensity distribution in the measuring area. 

The error analysis presented can be used to predict the accuracy in the calculated 
mean of the fluid dynamic quantities velocity, vorticity, turbulent intensity and 
Reynolds stress if the reading errors are known. 

Even using the present simple experimental set-up, and direct excitation of 
biacetyls phosphorescence at the border of its absoiption band (308nm), we can 
generate detectable emission with a observed lifetime of about 3ms. 
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Implementing phosphorescence of biacetyl vapor for the LIPA-technique requires 
a low oxygen concentration in the measuring area even if sensitizing of the 
emission by sulfur dioxide is used. A lower oxygen concentration than in the 
present experiments may improve the detectability of grid patterns. 

The use of sulfur dioxide to sensitize biacetyls phosphorescence is desirable if 
308nm excitation of the gas mixture is employed. The observed high absorption 
of incoming laser light by gaseous sulfur dioxide in the mixture requires more 
sophisticated concentration adjustment techniques. 
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Figure Bl: Taken pictures for set-up 1 


Figure B2: Taken pictures for set-up 2 
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Figure B3: Taken pictures for set-up 3 


Figure B4: Taken pictures for set-up 4 







Figure B5: Taken pictures for set-up 5 


Figure B6: Taken pictures without biaceyl 
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Appendix C: Data LIPA experiments 



Figure Cl: Numbering of the positions for mean 
turbulence intensity and Reynolds stress 
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Figure C2: Numbering of the positions for mean vorticity 
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2 

1 

0 
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-3 


distance from jet centerline, y/r distance from jet centerline, y/r 





nr 

position 

velocity [m/s] 


# 

<x/r> <y/r> 

<U X > <U V > <U> 



1.29 2.50 

-1.15 -3.32 3.93 

— 


2.37 2.51 

-0.49 -2.91 4.31 



3.58 2.38 

-0.71 -1.00 2.72 



4.58 2.49 

-1.99 -3.01 4.74 

^ 


5.56 2.28 

-1.66 -1.43 3.76 



6.66 2.28 

-1.68 -1.63 3.72 



7.10 1.41 

4.03 0.89 6.10 


8 

6.01 1.35 

4.80 0.26 5.99 


mm 

5.04 1.60 

3.96 0.43 5.03 


mm 

3.93 1.70 

1.40 0.31 3.28 


wM 

2.90 1.54 

1.40 0.18 3.35 


■9 

1.74 1.76 

! 1.63 -0.05 2.85 


mm 

0.65 1.79 

1.00 0.05 2.24 


mm 

1.18 1.03 

12.38 0.66 12.70 


15 

2.35 0.77 

14.78 0.41 14.94 


16 

3.34 0.86 

13.15 0.00 13.32 


17 

4.48 0.76 

13.20 0.03 13.45 


18 

5.55 0.67 

13.30 -0.08 13.63 

— 

19 

6.49 0.46 

12.48 0.28 12.73 


20 

6.90 -0.37 

12.38 0.08 12.62 


21 

5.99 -0.27 

14.68 -0.20 14.84 

i — 

22 

4.96 -0.16 

16.64 0.20 16.91 


23 

3.84 -0.11 

17.26 -0.10 17.40 


24 

2.73 0.07 

15.70 0.13 15.83 


25 

1.74 0.00 

16.16 0.00 16.32 


26 

0.55 0.26 

16.41 0.03 16.53 


27 

0.00 -0.39 

16.59 -0.26 16.61 


28 

1.15 -0.73 

15.21 -0.23 15.39 

— 

29 

2.13 -0.67 

14.98 -0.74 15.07 


30 

3.22 -0.87 

13.99 0.18 14.14 


31 

4.31 -1.04 

13.27 -0.10 13.44 


32 

5.40 -1.08 

12.15 0.03 12.36 


33 

6.39 -1.12 

11.59 -0.64 11.89 


34 

6.74 -1.91 

5.39 -0.56 7.34 


35 

5.78 -1.94 

4.54 -0.38 5.66 


36 

4.72 -1.98 

2.02 -0.84 4.09 


37 

3.67 -1.86 

1.30 -0.43 3.11 


38 

2.59 -1.64 

2.07 -0.23 3.57 

— 

39 

1.49 -1.44 

5.21 -0.77 5.80. 


40 

0.55 -1.42 

5.23 -0.87 6.15 


41 

0.92 -2.08 

0.79 -0.36 2.72 


42 

1.99 -2.40 

-0.79 0.05 2.34 


43 

3.11 -2.58 

-1.38 0.97 3.27 


44 

4.15 -2.73 

-1.43 1.89 3.64 


45 

5.15 -2.79 

-1.35 3.88 4.95 

- 

46 

6.16 -2.75 

-1.89 1.84 3.90 

L 

47 

1.53 -3.00 

1.30 3.06 4.34 


Table Cl: Velocities 
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CR [m/s] 

estimate [m/s] 




1.60 

0.25 

A 

-1.15 

3.71 

0.25 

A 

-0.49 

1.74 

0.25 

A 


2.59 

0.25 

A 

-1.99 

2.20 

0.25 

i A 

-1.66 

2.25 

0.25 

A 

-1.68 

4.53 

0.25 

M 

-1.61 

3.83 

0.25 



3.14 

0.25 

M 

2.57 

2.15 

0.25 

A 

1.40 

2.35 

0.25 

A 

1.40 

2.01 

0.25 

A 

1.63 

1.65 

0.25 

A 

1.00 

2.62 

0.25 

M 

3.21 

2.32 

0.25 

M 

-1.02 

3.26 

0.25 

M 

-0.64 

3.26 

0.25 

M 

-1.27 

3.77 

0.25 

M 

-1.58 

2.87 

0.25 

M 

-3.22 

1.99 

0.25 

M 

-3.48 

2.49 

0.25 

C 

-1.28 

2.58 

0.25 

C 

0.68 

3.12 

0.25 

c 

1.30 

1.92 

0.25 

c 

-0.26 

1.62 

0.25 

c 

0.20 

2.09 

0.25 

c 

0.45 

1.22 

0.25 

c 

0.63 

1.37 

0.25 

c 

-0.75 

2.33 

0.25 

c 

-0.98 

2.26 

0.25 

M 

0.38 

2.26 

0.25 

M 

3.18 

2.76 

0.25 

M 

2.61 

2.95 

0.25 

M 

2.42 

3.55 

0.25 

M 

4.79 

2.91 

0.25 

M 

4.52 

3.11 

0.25 

A 

2.02 

2.47 

0.25 

A 

1.30 

2.65 

0.25 

A 

2.07 

2.13 

0.25 

A 

5.21 

2.77 

0.25 

A 

5.23 

1.97 

0.25 

A 

0.79 

1.50 

0.25 

A 

-0.79 1 

2.75 

0.25 

A 

-1.38 

2.26 

0.25 

A 

-1.43 

3.07 

0.25 

A 

-1.35 

2.83 

0.25 

A 

-1.89 

2.93 

0.25 

A 

1.30 : 
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nr 

position 

turbulent 
intensity [%] 

1 

RCR [%] 

# 

<x/r> 

<y/r> 

u'/Uc 


(u'/U^RCR 

1 

1.29 

2.50 

9.44 

5.67 

14.17 

2 

2.37 

2.51 

21.85 

17.09 

6.53 

3 

3.58 

2.38 

10.23 

7.04 

13.45 

4 

4.58 

2.49 

15.26 

10.19 

8.90 

5 

5.56 

2.28 

12.93 

8.48 

10.56 

6 

6.66 

2.28 

13.21 

8.09 

10.36 

7 

7.10 

1.41 

26.62 

15.39 

5.25 

8 

6.01 

1.35 

22.55 

11.19 

6.10 

9 

5.04 

1.60 

18.46 

12.69 

7.43 

10 

3.93 

1.70 

12.63 

6.33 

10.91 

11 

2.90 

1.54 

13.82 

9.71 

10.05 

12 

1.74 

1.76 

11.80 

6.68 

11.77 

13 

0.65 

1.79 

9.69 

5.80 

14.40 

14 

1.18 

1.03 

15.40 

10.95 

8.56 

15 

2.35 

0.77 

13.66 

9.81 

9.62 

16 

3.34 

0.86 

19.18 

8.57 

6.88 

17 

4.48 

0.76 

19.15 

12.57 

6.89 

18 

5.55 

0.67 

22.20 

17.13 

5.96 

19 

6.49 

0.46 

16.90 

10.27 

7.80 

20 

6.90 

-0.37 

11.72 

7.68 

1131 

21 

5.99 

-0.27 

14.66 

9.50 

8.97 

22 

4.96 

-0.16 

15.17 

8.81 

8.66 

23 

3.84 

-0.11 

18.34 

10.44 

7.17 

24 

2.73 

0.07 

11.29 

6.94 

11.63 

25 

1.74 

0.00 

9.52 

637 

13.77 

26 

0.55 

0.26 

12.31 

7.82 

10.67 

27 

0.00 

-0.39 

7.20 

432 

18.21 

28 

1.15 

-0.73 

8.04 

530 

16.32 

29 

2.13 

-0.67 

13.68 

10.33 

9.60 

30 

3.22 

-0.87 

13.31 

933 

9.88 

31 

4.31 

-1.04 

13.27 

7.39 

9.91 

32 

5.40 

-1.08 

16.21 

10.79 

8.14 

33 

6.39 

-1.12 

17.37 

10.68 

7.60 

34 

6.74 

-1.91 

20.91 

12.74 

6.45 

35 

5.78 

-1.94 

17.12 

10.44 

7.90 

36 

4.72 

-1.98 

18.30 

12.34 

7.66 

37 

3.67 

-1.86 

14.50 

1035 

9.72 

38 

2.59 

-1.64 

15.60 

6.17 

8.96 

39 

1.49 

-1.44 

12.53 

835 

10.63 

40 

0.55 

-1.42 

16.32 

12.63 

8.23 

41 

0.92 

-2.08 

11.61 

734 

12.01 

42 

1.99 

-2.40 

8.80 

734 

15.64 

43 

3.11 

-2.58 

16.16 

10.82 

8.80 

44 

4.15 

-2.73 

13.31 

8.63 

10.31 

45 

5.15 

-2.79 

18.05 

14.87 

7.60 

46 

6.16 

-2.75 

16.64 

11.03 

8.38 

47 

1.53 

-3.00 

17.23 

11.43 

8.03 


Tabic C2: Turbulent intensities 
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position 


Reynolds stress 
[m 2 /s 2 ] 




RCR [%] 


< UV >RCR 


2.37 

19.48 

823 

12.60 

3.15 

16.29 

6.89 

12.34 

6.91 

12.95 

423 

14.02 

11.13 

9.52 

10.68 

9.84 

9.31 

10.52 

4.91 

13.52 

4.56 

11.95 

2.15 

16.06 

235 

19.11 

5.92 

11.56 

6.63 

13.84 

7.05 

11.93 

11.21 

10.73 

1222 

9.33 

1 6.94 

11.65 

5.13 

14.45 

5.90 

13.53 

8.13 

11.20 

829 

12.32 

425 

15.75 

4.03 

16.70 

4.67 

15.65 

1.43 

33.58 

421 

18.98 

2.83 

17.94 

4.01 

14.75 

5.39 

14.62 

5.77 

12.82 

10.21 

11.10 

36.42 

7.26 

10.60 

10.27 

13.91 

10.11 

533 

13.86 

8.64 

11.69 

5.68 

14.69 

1121 

11.22 

2.89 

15.95 

2.30 

18.25 

5.07 

14.00 

8.35 

12.85 

16.13 

12.07 

5.74 

13.15 

3.12 

17.38 


Table C3: Reynolds stress 
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nr position 


# 


1 1.22 1.77 

2 2.34 1.65 

3 3.44 1.62 

4 4.51 1.64 

5 5.54 1.47 

6 6.57 1.38 

7 6.01 0.56 

8 5.01 0.72 

9 3.90 0.80 

10 2.83 0.81 

11 1.75 0.89 

12 1.16 0.14 

13 2.24 0.04 

14 3.28 -0.01 

15 4.40 -0.14 

16 5.47 -0.21 

17 6.44 -0.32 

18 5.89 -1.10 

19 4.84 -1.06 

20 3.76 -0.97 

21 2.67 -0.78 

22 1.63 -0.71 

23 0.56 -0.57 

24 1.03 -1.42 

25 2.05 -1.54 

26 3.15 -1.74 

27 4.21 -1.90 

28 5.26 -1.95 

29 6.27 -1.93 

30 1.49 -2.23 


vorticity [10/s] 

[10* 1/s] 

<W Z > 

179.91 

72.09 

171.38 

79.96 

191.58 

75.93 

176.34 

84.52 

177.90 

92.78 

168.41 

84.03 

129.94 

90.68 

142.01 

82.18 

174.44 

71.99 

183.87 

64.53 

160.70 

62.89 

33.86 

49.68 

2.97 

62.94 

6.80 

54.47 

5.30 

51.14 

-24.61 

70.37 

-6.20 

69.82 

-134.25 

66.12 

-160.39 

72.88 

-185.82 

5458 

-140.78 

43.35 

-165.87 

69.68 

-129.34 

57.74 

-211.86 

96.60 

-176.43 

52.36 

-184.13 

61.34 

-180.98 

70.70 

-142.32 

94.64 

-167.58 

149.18 

-56.62 

66.63 



Table C4: Vorticity 
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Figure C4: Velocity contourplot 
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Figure Cl 1: Comparison of performed measurments and data from the literature. ▼: here 
performed measuments; A : measurments by Hussain and Clark [45] (Re d =360000); □: 
measurments by Crow and Champagne [46] (Re d =106000); x:measurments by Sami, 
Carmody and Rouse [47] (Re d =220000). 
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Figure Dl: Pictures taken with maximum intensification by the ICCD camera 
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Appendix E: Propert ies of chemicals 


r 

n 2 

S 02 

Biaceyl 

0 2 

1 M [g/mol] 

28.013 

64.06 

86.06 

31.999 

R [J/(kg*K)] 

297 

127 

95.6 

260 

p [kg/m 3 ] 

1.161 

(20° C; 1 atm) 

2.717 

(20° C; 1 atm) 

979.23 
(20° C; 
fluid state) 

1.326 

(20° C; 1 atm) 

| vapor 
flpressure [kPa] 



5.3 

(20° C) 


C [kJ/(kg*K)] 

I (20° C) 

1.038 

0.640 


0.917 

I ^ 

[kJ/(kg*K)] 
(20° C) 

0.741 

0.504 


0.656 

K 

1.401 

1.27 


1.398 I 


Table El: Properties of chemicals [15,49,50,51] 
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required statistical safety [%] 


■SH 

95.0% 

99.0% 

99.9% 

i 

12.71 

63.66 

636.62 

2 

4.30 

9.93 

31.60 

3 

3.18 

5.84 

12.94 

4 

2.78 

4.60 

8.61 

5 

2.57 

4.03 

6.86 

6 

2.45 

3.71 

5.96 

7 

2.37 

3.50 

5.41 

8 

2.31 

3.36 

5.04 

9 

2.26 

3.25 

4.78 

10 

2.23 

3.17 

4.59 

20 

2.09 

2.85 

3.85 

30 

2.04 

2.75 

3.65 

40 

2.03 

2.71 

3.55 

50 

2.01 

2.68 

3.50 

60 

2.00 

2.66 

3.46 

70 

1.99 

2.65 

3.44 

80 

1.99 

2.64 

3.42 

90 

1.99 

2.63 

3.40 

100 

1.98 

2.63 

3.39 

200 

1.97 

2.60 

3.34 

300 

1.97 

2.59 

3.32 

400 

1.97 

2.59 

3.32 

500 

1.97 

2.59 

3.31 

1000 

1.96 

2.58 

3.30 

oo 

1.96 

2.58 

3.29 


an as a function of number of 
required statistical safety [ 16 ]. 








